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Computer Software Management 
and Information Center 


The documentation and program developed for the 
advanced structural geometry studies will be 
made available to the public through COSMIC. 

COSMIC (Computer Software Management and Information 
Center) was established early in 1966 at the 
University of Georgia to collect and disseminate 
to the public computer software developed by govern- 
ment agencies. Since that time thousands of computer 
programs in all areas of aerospace engineering, mathe- 
matics, business, and industry have been distributed 
to requesters throughout the United States. 

The Technology Utilization Division of NASA, 
designed to enlarge the return on the public invest- 
ment in aeronautical and space activities, was the 
first government agency to participate formally. 

In July 1968 the Atomic Energy Commission and in 
November I968 the Department of Defense Joined in 
the COSMIC endeavor. With the addition of these two 
major agencies, the original concept of making tax- 
paid developments available to the public was expanded 
to make COSMIC a transfer point between and within 
government agencies as well. 

Requests for documentation or information concerning 
this program should be directed to: 

COSMIC 

The University of Georgia 

Barrow Hall 

Athens, Georgia 306OI 

REF: HQN- 10677 ' 
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1.1 INTRODUCTION 


One of the most economical structural systems in 
contemporary use has been based on the spherical form. 
Designs for such structures are influenced primarily 
by the ultimate purpose of the structure, special 
environment to which it will be subjected, and the 
materials of fabrication. 

Two basic systems are used for subdividing the 
spherical form for structural application: The bi-polar 

system and the multi-polar system. Figures 1.1, 1.2. 




I 

I 


Bi-polar System 


Multi-polar System 


Figure 1.1 


Figure 1.2 
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The bi-po1ar system 1s related to the familiar 
latitude-longitude approach to subdividing a sphere. Two 
common examples of this system are the ribbed dome (Figure 
1.3) and the lamella dome (Figure 1.4). 




Ribbed Dome 
Figure 1.3 
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Lamel 1 a Dome 


The multi-polar system is related to the spherical form 
of polyhedra. Perhaps the most familiar example of this 
system is the geodesic dome discoverd by R. Buckminster 
Fuller. Figure 1.5 



Geodesic Dome 
Figure 1.5 


A typical design problem with the bi-polar system is 
with the geometrical relationships of the frame configura- 
tion. Computer aids in the design investigation may be 
initiated to handle the great number of variables in the 
determination of lengths, frequency of members, total 
number of joints, number of members intersecting at a 
single joint, relationship of members to each other, etc. 

This portion of the report will concern itself with 
a mathematical and computer model for determining the 
geometrical properties of the multi-polar system for sub- 
dividing a spherical form for structural applications. The 
model has been limited to the three polyhedral forms made 
up completely of regular triangles: the tetrachedron , 

octahedron and icosahedron. 
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1.2 POLYHEDRON 


Hoppe, in 1882, coined the word polytope: a geo- 

metrical figure bounded by portions of lines, planes, or 
hyperplanes; in two dimensions it is a polygon; in three a 
polyhedron.* However, the Greeks studied polyhedra over 
two thousand years ago with the findings of Euclid. Others 
such as Klein, Schiafli and Coxeter, introduced much to the 
study and concepts of the polytope. In this section of 
the report, the structural configurations discussed are 
based on several of the polyhedral forms, specifically the 
regular (Platonic) polyhedra; the Tetrahedron, Octahedron, 
and Icosahedron. 

In Euclid's writings. The Elements , explanation and 
definition is give to the five regular solids as known 
to the ancient world. The convex polyhedra are said to 
be regular if each have regular and equal faces, if they 
are congruent, and if they are of regular polyhedral angles. 
Table 1.1 list the properties of the Tetrahedron, Octahedron, 
and Icosahedron which are considered as three of the five 
regular polyhedral forms. 


*Coxeter, N. S. M. 1 



Tabl e 1.1 


PropsrtiGS of the Basic Polyhedra 
Tetrahedron, Octahedron, Icosahedron 


Te trah edro n 


3 


3 


V=4, F=4, E=6 


Dihedral angle/3= 2 sin/y/3 - 70° 31' 44 


Angle subtended by an 
edge at center of 


polyhedron 


Verti ces 


Center to mid edge 
center to center of face 
center to vertex 
edge 

mid edge to center of face 

mid edge to vertex 

mid edge to opposite mid edge 

height (vertex to center of 
opposite face) 


5 = cos(- 1/3) = 109°28'16" 

(1/yr, 1//T, 1/yT) (i//r, -i//r, -v/ 3 ) 

(-i//r, i/yr, -i/tt) (-i//y, v/^) 

i/yr 


1/3 

1 

isTEiJz 
727 3 
/2 

4/3 


0.57735 

0.33333 

1 

1 .63299 
0.471 40 
1 .41421 
1 .1 5470 
1 .33333 


area of face 
volume 


2//y 1.1 5470 

8/y 0.51320 

27 
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Table 1.1 (cont.) 


Octahedron 3 
V=6, F=8, E=12 

Dihedral Angle)3= tan Ijl = 109° 28' 16" 

Angle substended by 
an edge at center 
of polyhedron 8= 90° 

Vertices (±1 , q, 0) 

(0. il , 0) 

(0, 0, il) 


edge 


1 .41 421 

center to vertex 

1 

1 

center to mid edge 

1//2 

0 . 70711 

center to center of 
face 

i/yr 

0 .57735 

Mid edge to near 
vertex 

/ 3/2 

1 .22474 

Mid edge to distant 
vertex 

/ 5/2 

1.58114 

area of face 

/2 

0.86603 

Vol ume 

4/3 

1 . 33333 
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Table 1 . 1 (cent . ) 

5 

Icosahedron 3 


V=12, F=20, E=30 

Dihedral angle e= t -(1 sin ) = 138°11'22 

Angle subtended by 

an edge at center of _ ^ oofrinc 

polyhedron 6= cos (/5/5 ) - 63 26 05.81 8 


Vertices t - 1 - 1.61803 

2 


("■ 


+ \/^ 


1 


1/4 

5 W T 


1 = T - 1 

/+ 1 

, 0 , + '/T 

X 

V bV-'/t 


= X + 1 



x/f = X + 2 

A 7T 

+ 1 

V 5^/^ 

5'''VX 

edge 

2/5^^'' /T 

1.05146 

center to vertex 

1 

1 

center to mid edge 

1/4 

/7/5 ' 

0.85065 

center to center of 
face 


0 . 79465 

area of face 

1/ts/^ 

0.47873 

vol ume 

(4(5 / ^ ) y~r ) / 3 

2.5361 5 
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1.3 STRUCTURAL ORIENTATION 


The structural configuration desired is acquired 
through a three-way gridding of the faces of the polyhedral 
form chosen from one of the three triangular faced regular 
polyhedra. The grid is then translated onto the surface 
of a circumscribed sphere. A three dimensional rectangular 
coordinate system was chosen for the basis of the compu- 
tations. Due to the symmetries existing in the polyhedral 
forms only one face of the polyhedron is used in calculation 
of the geometrical properties of the structure. Figure 
1.6 shows the orientation of the polyhedral form with 
respect to the x, y, z axis with the vertices of the face 
chosen for the geometrical computations. The intersection 
of the X, y, z axis is located at the origin (0,0,0) of 
the polyhedron, this point being the center of the cir- 
cumscribed sphere. Table 1.2 list the coordinates of the 
vertices of the faces chosen as the PPT. 



Octahedron Tetrahedron 

Figure 1.6 
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Icosahedron 
Figure 1.6 ( con t . ) 

Table 1.2 

Coordinates of the Principal Polyhedral Triangles 

T etrahedron 

= (-1//T, -1//3", 1/yi) = (-.57735027 , -.57735027 . .57735027 ) 
P 2 = (1//^. -V/T, -1 //3) = (.57735027 . -.57735027 , -.57735027) 
P3 = (-1//T, l/A, -1 /-/3) = (-.57735027 , . 577 35027 .-. 57735027 ) 

Octahedron 

Pi =- n , 0 , 0 ) 

P 2 = ( 0 , 1 , 0 ) 

Pj, = ( 0 , 0 , 1 ) 

I cosahedron 

Pi = 1/5 / 7 T) = (0 , .85065081 , .525731 1 1 ) 

P 2 = ( 1 / 5 ^/VT 0, /T/5^/'') - (.525731 1 1 , 0 , .85065081 ) 

P 3 = (x/7/5^/'', 1/5^/''/7", 0 ) = (.85065081 , .525731 1 1 , 0 ) 
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Throughout the discussion of the methods of sub- 
division of the polyhedral forms the Icosahedron will be 
used with examples of computer maps of the spherical forms 
derived using the three traditional orientations: edge, 

vertex, and face. Figure 1.7. 



Edge Orientation 
Figure 1.7 
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Vertex Orientation 

Figure 1.7 (cont.) 


1 - 11 + 



-L 


Face Orientation 

Figure 1.7 
( con t . ) 
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1 .4 DEFINITIONS 


AXIAL ANGLE (fj) = an angle formed by an element and a 

radius from the center of the polyhedron meeting in a 
common point and the vertex of the axial angle sharing 
a vertex of the polyhedron. 

CENTRAL ANGLE ( 6 ) = an angle formed by two radii of the 

polyhedron passing through the end points of a principal 

side. 

CHORD FACTOR (cf) = the element lengths calculated based 
upon a radius of a non-dimensional unit 1 for the 
spherical form. The length of any element for larger 
structures may be found by the equation; 

cf X r = 1 

where cf = chord facto>^ 

r = the radius of the desired 
structural forms 

1 = the length of the element 
sought 

DIHEDRAL ANGLE ( 3 ) = an angle formed by two planes meeting 
in a common line. The planes themselves are the faces 
of the dihedral angle, and the common line is the element. 
To measure the dihedral angle measure the angle whose 
vertex, is on the element of the dihedral angle and whose 
sides are perpendicular to the element and lie one in 
each face of the dihedral angle. 
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FACE ANGLE (a) - an angle formed by two elements meeting 
in a common point and lying in a plane that is one of 
the faces of the polyhedron. 

FACES = the triangles making up the "exploded" structural 
form . 

FREQUENCY (v) = the number of parts or segments into which 
a principle side is subdivided. 

PRINCIPLE POLYHEDRAL TRIANGLE (PPT) = any one of the equal 
equilaterial triangles which forms the face of the 
regular Polyhedron. 

PRINCIPLE SIDE (PS) = any one of the three sides of the 
principle polyhedral triangle. 



1 .5 METHODS OF SUBDIVISION 

Upon using the spherical form as a structural unit, 
it is readily apparent that the basic polyhedral form, 
in its pure state, can not satisfy the range of conditions 
that must be geometrically and structurally met. Seven 
methods will be discussed for reducing the basic poly- 
hedral form into a larger number of components from which 
the geometrical properties may be made to remain within the 
structural fabrication and erection limits for a desired 
conf i gurati on . 

Due to the symmetrical characteristics of the basic 
polyhedral form only one face of the polyhedron is used 
for calculating the geometrical properties of the struc- 
tural configuration. The remaining faces may be found by 
rotations or reflections of the principal polyhedral tri- 
angle and its transformations. 

Attention is given here to the seven methods of sub- 
dividinq the PPT in a broad sense and will be treated in 
detail in the following section. 

Method 1 : 

The PPT is subdivided into n frequency, with the parts 
chos.en as equal divisions along the three principal, sides. 

Fi gure 1 . 8 ^ 

NOTE: FT = T2 


Figure 1.8 
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Each point of subdivision is then connected with a 
line segment parallel to their respective sides thereby 
giving a three-way grid so that a series of equilateral 
triangles are formed. Figure 1.9 



Note: AB is parallel to 12 


Each vertex on the PPT is then translated along a line 
passing through the origin (0,0,0) of the polyhedron and 
its respective vertex, onto the surface of the circumscribed 
sphere. The element connecting the translated vertices form 
the chords of a three-way great circular grid. Figure 1.10 



Figure 1.10 
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Methods 2 & 3: 


The PPT is subdivided into n frequency with the parts 
therein as equal arc divisions of the central angles of 
the polyhedron. Figure 1.11. 


A 



Note; A1 f 12 


The points of subdivision on each principal side of 
the PPT are connected with line segments parallel to their 
respective sides. Each line segment intersects at a number 
of points which define a grid of subdivision. Due to the 
method of subdivision, small equilateral triangular "windows 
occur in the grid. Figure 1.12. 

1 A 



NOTE : AB i_s paral 1 el 

to 12 


A a ^ a b 

Windows are equilateral 
tri angles 
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The center of these "windows" are found by one of two 
methods and are used as the vertices of the three-way grid 
for the PPT. They are then translated onto the surface of 
the circumscribed sphere along a line passing through the 
respective vertex and the origin (0,0,0) of the polyhedron. 
The element connecting the translated vertices form the 
chords of a three-way great circular grid. Figure 1.13. 



Figure 1.13 
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Method 4: 


The PPT is subdivided into n frequency, with the parts 
chosen as equal divisions along the three principal sides. 

Fi qure 1 . 1 4 . ^ 

1 

2 

Note: A1 = 12 


Figure 1.14 

Each point of subdivisions is then connected with line 
segments perpendicular to their respective principal side 
thus giving a three-way grid comprised of equilateral and 
right triangles. Figure 1 .15. 





Figure 1.15 


Each vertex on the PPT is then translated onto the 
surface of the ci rcumscri bed sphere along a line passing 
through the respective vertex and the origin (0,0,0) of 
the polyhedron. The elements connecting the translated 
vertex form the chords of a three-way great circular grid. 
Figure 1.16 



Figure 1.16 
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Method 5 

The PPT is subdivided into n frequency with the parts 
chosen as equal arc divisions of the central angle of the 
polyhedron. Figure 1.17. 

NOTE: T\ ^ T2 


Figure 1.17 

The points of subdivision on each principal side of the 
PPT are connected with line segments similar to Method 4. 
However, the line segments are not perpendicular to their 
respective sides. Upon completion of the connections a 
grid is created. Due to the method of subdivision, small 
triangular "windows" occur in the grid. Figure 1.18. 

NOTE: M 

Small triangular 
windows occur 

Figure 1.18 
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The centers of these "windows" are found and are used 
as the vertices of a three-way grid for the PPT. The ver- 
tices are then translated onto the surface of the circum- 
scribed sphere along a line passing through the respective 
vertex and the origin (0,0,0) of the polyhedron. The 
elements joining the translated vertices form the chords 
Of a three-way great circle grid. Figure 1.19. 



Figure 1.19 
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Method 6: 


The PPT may be described as six riqht triangles each 
being a reflection or rotation of the other. Figure 1.20. 

Note: ABC is 

a right 
tri angl e 


Figure 1.20 

In this method of subdivision we shall treat only triangle 
ABC. The remaining section of the PPT may be found through 
rotations and reflections of this basic unit. The Line AB 
is subdivided into parts chosen as equal arc divisions of 
the central angle of the polyhedron. Figure 1.21. 

Note: A1 f 12 



A 
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Once the subdivisions are found they are used to find 
the points of division on side AC and CW. Perpendiculars 
through the points of division on side ^ are extended to 
side AC, this giving the points of subdivision on side 
AC . Figure 1.22. 


Note: 24 J_ AB 
AT ^ 23 


Figure 1.22 

The points of division on the side CB" were formed by 
extending a line through the points of subdivision on 
side AC perpi ndi cul ar to side Figure 1.23. 


c 



Note : 54 _L CB 

T2 ^ JA 

c 



Figure 1,23 
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Having acquired the points of subdivision along the 
three sides of the triangle, diagonals are drawn from each 
point on side AC to alternate points of sides AB and BC . 
Figure 1.24. 


c 



Fi q ure 1 . 24 

To complete the three-way grid connect alternate 
points of subdivision of side AB to alternate points of 
division of side FC. Figure 1.25 


c 



Figure 1.25 



Through rotations and reflections of the basic unit 
and its subdivisions, the entire three-way gridding of the 
PPT may be found. Figure 1.26 



Figure 1.26 

The vertices of the three-way grid are then trans- 
lated to the surface of the circumscribed sphere along a 
line passing through the respective vertex and the origin 
(0,0,0) of the polyhedron. The element joining the trans- 
lated vertices form the chords of a three-way great circle 
grid. Figure 1.27. 
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Figure 1.27 
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Method 7 


The PPT is described as six right triangles each 
being a reflection or rotation of the other. Figure 1.28. 



NOTE: ABC is a 

right 
triangle 


Figure 1.28 

In this method of subdivision only one of the right 
triangles will be treated as the basic unit for subdivision 
into a three-way grid. The line is subdivided into parts 
chosen as equal arc division of an angle made up of the triangle 
AC and the origin of the polyhedron with the origin (0,0,0) 
bieng the center of the triangle of subdivision. Figure 1.29 



NOTE: A1 ^ 23 

WC ±IC 


Figure 1.29 
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Once the subdivisions are found on line AC they are 
used to find the points of division on side AB and CB. The 
lines through the points of sibudvision on side AC are taken 
perpendicular to side A¥, this giving the points of division 
on side A¥. Figure 1.30. 



Figure 1.30 

The points of division on CB are found by extending a 
line through the points of subdivision on XC perpendicular 
to ^B". Fi gure 1.31. 



Note: 54 _L CB 

T2 5T 


Figure 1.31 



Having acquired the points of subdivision along the 
three sides of the triangle, diagonals are drawn from each 
point on AC to alternate points on AB and BC. Figure 1.32. 



Figure 1.32 

To complete the three-way grid connect alternate points 
on M" to alternate points on bT. Figure 1.33. 


c 



Figure 1.33 


Through rotations and reflections of the basic unit 
found and it's subdivison, the entire three-way grid of 
the PPT may be found. Figure 1.34. 



F i gure 1.34 

The vertices of the three-way grid are then translated 
to the surface of the circumscribed sphere along a line 
passing through the respective vertex and the origin (0,0,0) 
of the polyhedron. The elements joining the translated ver- 
tices form the chords of a three-way great circle grid. 
Figure 1.35. 


1 - 5 ^ 



Figure 1.35 
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1 .6 METHOD 1 


The mathematical and computer model was developed 
for subdividing a tetrahedron, octahedron, or an icosahedron 
circumscribed by a unit sphere. The icosahedron was chosen 
as an example to illustrate the geometry of the model. 

The icosahedron is oriented in a three dimensional 
rectangular coordinate- system so that the vertices of one 


P P T are: 

( X 1 , y I i Z I ) 


(Xj, Yj, Z^) 


= ( 0 a . 

Vi 

= (0, .850651, 

= f 1 0, 

^ ( .525731 , 0 , 


__J ) 

4 

nAn/V 

.525731 ) 



.850651 ) 


(X3, Y3, Z3) 


( /T , 1 0) 

4 4 

y ^/7 

( .850651 , .525731 , 0 ) 


where x 


= 1 +A" 


2 

with the intersection of the axis X, Y, Z located at the 
origin (0,0,0) of the icosahedron. Figure 1.36. 
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{ 0,0 ) 


Figure 1.37 
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To find the projection of each vertex of the PPT 
onto the unit sphere along a line segment through the vertex 
of the PPT and the origin each coordinate of each vertex, 
PPT, is divided by the distance between the vertex PPT and 
the ori gi n . Figure 1.38. 



Figure 1.38 


Using the coordinates, this program finds the lengths 
of the elements of the structure («.), the angle between pairs 
of elements (face angle c), the angle between the elements 
and a radius from the origin to an endpoint of the element 
(axial angle n) , and the angle between adjacent faces of the 
structure (dihedral angle b)- Figure 1.39 



Fi gure 1 . 39 

To find the angle between elements the face ) a, we use 
the coordinates of their endpoints. The vertex of the angle 
is a common endpoint to each element and is translated to the 
origin. The other two endpoints Pi and ?2 are translated in 
the same manner. Letting (Xi, Yi, Zi) and (X 2 , Y 2 , Z 2 ) be 
the points resulting from the translations of the endpoints 

?i and P 2 9 


I-ilO 



cos a 


XjX, + YjY^ + ZjZ, 


where di 
and d2 

a is the desired angle. 

To find axial angles the above method is used except 
that the vertex is established at one end of an element and the 
origin is used with the other endpoint to define the angle. 

The desired angle is n. 

The angle between two adjacent faces, the dihedral 4 3, 
is found using 

cos 3 = AjA2 + B1B2 + C1C2 

222222 

Ai + Bi + Cl Aj + + Cj, 1.3 

where 

3 is the desired angle. 

AiX + BjY + CjZ + = 0 defines the plane containing 

one face and A2X + B2Y + C2Z + D2 = 0 defines the plane con- 
taining the other face. 

The negative sign is used because the obtuse angle is desired. 
The A, B, and C for each plane are computed as 
Vi Zi 1 

A = Y2 Z2 1 
Y 3 Z3 1 
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B 



Xi 

Zi 

1 

B = 

X2 


1 


^ 3 


1 


Xi 

Yi 

1 

C = 

X2 

Y2 

1 


X3 

Y3 

1 


1 .4 


where (Xi, Y^, Zi), (X2> Y^» ^2^’ (Xs» Y3, Z3) lie in 

the plane. In particular the three vertices of each face 
are used. For the special case where the two faces used lie 
above separate faces of the polyhedron, the assumption is 
made that the plane containing the element common to each 
face and the origin bisects the angle. This angle is found 
in the same manner and doubled. This method is used because 
the structural face above the polyhedral face is not generated 


properly . 

The length of the elements, , are found by using the 



£ is the desired length. 

To reduce total output, this program takes into account 
certain symmetries and outputs only a part of the total angles 
and lengths. The rest of the values are the same as at least 
one outputed value and can easily be found using the following 
symmetries. Figure 1.40. 
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F1 gure 1.40 
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FACE ANGLES 

For every face angle opening directly towards (or 
away from) the point (0.0), there are egual angles opening 
towards (or away from) the point (N.O) and (N.N). For 
example, the angle (1.1). (0.0), (1,0) with vertex at (0.0) 
is egual to the angle (N-1,0). (N.O), (N.l) and the angle 
(N,N-1), (N.N), (N-1, N-1). Thus, only the face angles 
facing directly towards or away from (0,0) are computed. 

If the vertex is to lie at (I.J), the angle will be either 

(1 + 1,0 + 1), (1,0), (I +1. 0) (I-l, 0-1). (I.O). 

(1-1,0). Also, only the face angles falling on the right 
of or on a line passing through (Xi, Yi, Zi,) and the mid- 
point of the opposite side are computed. 

The elements of the structure can be put into one-to- 
one correspondence with the lengths and dihedral angles. 

The dihedral angle associated with an element is the angle 
between the two faces containing the element. For each 
element, there are two axial angles, one at each end, but 
since the element is a cord of the circle, the two angles 
are equal and may be considered one. In this case, we 
have a one-to-one correspondence between elements and axial 
angles. This programwill only compute values around 
elements parallel to the side opposite (Xj, Y^, Z^) and on 
the right side of a line through (X^, Y^, Z^) and the 
midpoint of the opposite side. All other lengths and angles 
are symmetric to one of the lengths and angles computed in 

this manner. 
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THE COMPUTER PROGRAM DESCRIBED ON 
PAGES 1-45 to 1-75 


IS AVAILABLE FROM COSMIC 



1 .7 METHODS 2 & 3 


This program works with a tetrahedron, octahedron, 
or icosahedron circumscribed by a unit sphere. The icosa- 
hedron was chosen as an example to illustrate the geometry 
of the program. The icosahedron is oriented in a three 
dimensional rectangular coordinate system so that the 


vertices of one PPT are 


X 
1— ■ 

I—* 

N 
1— ■ 

M 

( 0 , , 1 ) 
V — r 

\/5 v5\/t 


(0, .850651, .525731) 

^ ^ 2 » 2 ’ ^ 2 ^ 

(1 , 0, /r ) 


( .525731 , 0 , .8506 51 ) 

( X 3 , y 3 , z 2 ) 

( \/t , 1 » 0 ) 

\/5 V 5 V T 


(850651 , . 525731 , 0 ) 

where x " 

1 + \/^ 

2 

with the intersections of 

the axis X, Y, Z, located at 


the origin (0,0,0) of the icosahedron. Figure 1.43. 





z 



This PPT is divided into smaller triangluar units which 
are translated onto the surface of a sphere constituting 
the desirable space form. 

Using the following formula the planes consisting of 
the edges of the PPT and the origin (Xj, Y^, Zj) (X2, ^ 2 ^ ^2^ 
(X3, Y3, Z3) are rotated from 3 -space into 2 -space, Figure 
1 . 44 . 


I-TT 









Figure 1.44 

x' = XiX + Uj^y + v^z 1.6 

y' = x.x + y y + V z 
^22 

z ' = A3X + + V3Z 

Where X, y, v are direction cosines of the X'-axis, and Y'-axis, 
and Z'-axis respectively with respect to the old axis and 
are found by: 



Xzs X 3 ; y 2 »U 3 ; andv 2 ,v 3 are found similarly. 
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cu 


(U 


y 1 ( X 1 - X 2 ) + X 1 ( y 2 - y 1 ) 


OP i s 


o 



1-80 



C2 

D 

Rotate the points of intersection along the PPT edge from 
2 -spaces back to 3 -spaces. 


X 

II 

X 

+ Mi-y' 

+ v,z 


1.11 

y = A 2 X' 

+ May' 

+ V2^ 



X 

CO 

II 

N 

* ysy' 




wherex, y» v. dire 

direction cosines 

of the X' 

'-axis, Y'-axis, 

and Z'-axis with 

respect 

to the old 

axis and 

are found: 

f— 1 

X 

It 

rH 

<< 

\/xi^ + 

/ J‘ 

2 ^ 2 

Vl + Zj 



VI ^ yi / 


2 2 
.Vi + Zj2 



''1= z, / 

y>^ 

2 ^ 2 
y 1 + z 1 




^2» ^ 3 » 1^2* 1^3^ ^2» ^3 found Similarly. 

Retain the co-ordinates along the edges S^, $2 and S3 as 
shown i n Fi gure 1.47. 
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Figure 1.47 

After finding the unit measurements along the edges of 
the PPT, they are connected through a grid determining a 
smaller grid network. Since the units along the PPT edge 
are not equal, the gridding will create "windows". The 
centers of these "windows" must be found to establish the 
final 3-way grid network. Figure 1.48. 



Figure 1.48 
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The gridding and windows are found by the following 
method: From the coordinates along the edges of S^, S2 

and S3, calculate coordinates of the window by finding the 
intersection of P2 with Ps" >4 and PT~^ with 
P5 Pe and P3 P4 with P5 Pe using the two point form of 
the eguation of a line in three-space for the three lines 
and solve simultaneously for the points of intersection. 


P1P2 

i s : 

X - 

Xl 

y 

“ yi 

^ 



X 2 - 

Xl 

^2 

- 

N 

1 

CM 

N 

*^3^4 

i s : 

X - 

X 3 

y 

- V3 

Z - Z 3 



X 4 - 

X 3 

^4 

- ^3 

z - z 

4 3 

^5^6 

i s : 

X - 

X 5 

y 

- >5 

Z - Z g 



Xg-X 

5 

^6 

-^5 

^6- ^5 


To find the intersection of P^P^ with PgP^^ the equation 
takes the following form: 

(1) P1P2 is: x(y2-yi) + y(xi-X2 ) = yi(xi-X2) + XiCy^-y^) 

(2) is: y(z^-Zj) + z(,V2-yi ) = 

( 3 ) FTFi, is: x(y4-ys) + y(xs-X4 ) = yi(x3-X4) + Xj(y4-y3) 

( 4 ) P7P4 is: y(z4-Zj) + z(y3-y4 ) = Z3(y3-y4) + ^3(24-23) 
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For P1P2 (y2~^i) 

(Xj-Xj) = bj 

yi{Xi-x^) + Xj(y2-yj) = Cj 

For P3P4 let: ^ ®2 

(Xj-x^) = 

using the formula 1.10 solve for x and y coordinates of 
the intersections of Pi^a with P3P4- 
Find the z coordinate: 

For P1P2 

(Yi-Ya) = 

Ziiyi-Vz'^ yi(z2"23) " 


For P7P4 (24-23) 

(Y 3 -y 4 ) " 
Z3(y3-Y4) 

The other two vertices of the 
manner. Once the coordinates 
are determined, its center is 


»2 

^^2 

+ y3(24’*Z3) = ^2 

window are found in a similar 
for the vertices of the window 
found by one of the following 


two methods : 

METHOD I: 

On the PPT Plane the windows appear as equilateral 
triangles with vertices Pi(XiyiZ^), P2^^2^2^2^’ 
P‘(x3y3Z3) as shown in Figure 1 . 49 . 
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Figure 1.49 


The center C(cx,cy,cz) is found with the following 
formul a : 

CX = Xi+ X 2 + X 3 1.13 

3 

CY = .Y1+ y^+ Y3 


cz 


3 


^1+ 2 2+ Z3 


METHOD II: 


3 


The coordinates of the window found on the surface 
of the PPT are first "exploded" to the surface of 
the sphere. The center of the exploded window is 
then found by the intersection of angle bisectors. 
To find the projection of each vertex of the window 





onto the unit sphere, translate each vertex along 
a line through the vertex of the PpT and the origin; 
each coordinate of each vertex, PPT, is divided by 
the distance between the vertex PPT and the origin, 
Figurel.50. 



Figure 1.50 
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2 



where d = distance 1.14 
from origin to P 

1 

where r = radius of 
the sphere to be ex- 
ploded upon and is 
considered constant 


X ' 1 = rx 

1 

^ 1.15 

y'l = ry^ 

“d~ 



z 


Figure 1.51 
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PiTx = Pi. - 


P^Ty ” P'y ' P^y 
P^Tz = P ^ - 


P^Tx 


= Pz - Ps„ 

X X 


p - = P 2 “P 3 

^2Ty X 

P? =P9 - P 3 

^^Txi ’ ^^Tyi ’ ^^Tzi 


= 0 


Rotate plane Pi, P 2 . ^3 so that P 1 P 3 will fall on 
and P 3 is at the origin using equation 1.6 

The center is found with the intersection of 
bisectors of the triangular window PiPaPs* Figure 



Figure 1.52 


1 .16 


the X-axis 

two angle 
1.52. 
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The angles y and 5 are found: 


Arctan 1.17 

= Y 

Arctan 

= 6 

Xi-X^ 

rotate ?2 about Pg toward P^, ^/2 y degrees 

x^ = X2 cos V2 Y + Y2 sin 1/2 Y 1.18 

y4 = y2Cos 1/2 Y - x^ sin V2 y 

locate the origin at P^, then rotate P^ about Pj toward 
Pg ^ / 2 ^ degrees. 

X5 = (X2-X1) COSV2 < 5 - y^ sin V2 <5 + Xj 1.19 

y 5 ~ y 2 cos ^ / 2 5 + (x2“ X;^) sin V2 ^ 
thus defining Pi P5 and P3P4 . 

With P3 at the origin formula 1.9 mav be used to solve 
for the intersection of line ^^4 , P7P5 finding center C. 
Rotate C back to three space using formula 1.11 . Then 

translate center C back to three space ("C" is located in 
the previously "exploded" window), figure 1.53. 

C ' X = Cx + P3X 1.20 

C'y = Cy + P3.y 
C'z = Cz + P3Z 
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Fi q u re 1.53 

For Method I or Method II, the centers found are 
"exploded" to the surface of the sphere using formula 
1.14 and f ormul a 1.15 

Using the coordinates, this program finds the lengths 
of the elements of the structure {a), the angle between 
pairs of elements (face angle a), the angle between the 
elements and a radius from the origin to an endpoint of 
the element (axial angle fi), and the angle between adjacent 
faces of the structure (dihedral angle g), Figure 1.54. 
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Fi qure 1.54 


To find the anqle between elements the face ^ a, we 
use the coordinates of their endpoints. The vertex of the 
anqle is a common endpoint to each element and is translated 
to the origin. The other two endpoints Pj and P 2 are trans- 
lated in the same manner. Letting (xj, yj , zj and (xj, yz, 
Zj) be the points resulting from the translations of the 
endpoints P^ and P 2 , 


cos a 

where 
and d^ 
a is the desired angle. 

To find axial angles the above method is used except 
that the vertex is established at one end of an element 
and the origin is used with the other endpoint to define 
the angle. The desired angle is 

The angle between two adjacent faces, the dihedral 
^ 3 , is found using 




where 

3 is the desired angle. 

+ B^Y + C^Z + Di = 0 defines the pi ane conta i ni ng 
one face and A^X + B^Y + C^Z + D 2 = 0 defines the 
plane containing the other face. The negative sign is 
used because the obtuse angle is desired. 

The A,B, and C for each plane are computed as 

A = Y Z 1 

2 2 

3 ^3 ^ 


1 .23 
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Xi Yi 1 

C = 

X 2 ^2 1 

X 3 Y 3 1 

where (xi, Y , Z^), (X 25 Y 2 » Z 2 )» and (X 3 , Y 3 , Z 3 ) lie in 
the plane. In particular the three vertices of each face 

are used. 

The length of the elements^ are found by using the 

general 
i = 

z is the desired length. 

To reduce total output, this program takes into 
account certain symmetries and outputs only a part of the 
total angles and lengths. The rest of the values are the 
same as at least one outputed value and can easily be 
found using the following symmetries, Figure 1.55. 


egua ti on : 


(P - P )^ + (P - P )^ + (P^ - Pf 


1 .24 
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FACE ANGLES 

For every face anqle opening directly towards (or 
away from) the point (0,0), there are equal angles opening 
towards (or away from) the point (N,0) and (N,N). For 
example, the angle (1,1), (0,0), (1,0) with vertex at (0,0) 
is equal to the angle (N-1,0), (N,0), (N,l) and the angle 
(N,N-1), (N,N), (N-1, N-1). Thus, only the face angles 
facing directly towards or away from (0,0) are computed. 

If the vertex is to lie at (I,J), the angle will be either 
(I + 1, J + 1), (I,J), (I + 1, J) or (I-l, J-l), (I,J), 

(I-1,J). Also, only the face angles falling on the right 
of or on a line passing through (X^, Yj, Z^,) and the mid- 
point of the opposite side are computed. 

The elements of the structure can be put into one- 
to-one correspondence with the lengths and dihedral angles. 

The dihedral angle associated with an element is the angle 
between the two faces containing the element. For each 
element, there are two axial angles, one at each end, but 
since the element is a cord of the circle, the two angles 
are equal and may be considered one. In this case, we have 
a one-to-one correspondence between elements and axial 
angles. This program will only compute values around elements 
parallel to the side opposite (X^, , Zj and on the right 

side of a line through (X^, Y^,Zj) and the midpoint of the 
opposite side. All other lengths and angles are symme t r i c 
to one of the lengths and angles computed in this manner. 
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The computer program here contained was written 
for the IBM 7040/7044 computer, utilizing FORTRAN IV 
language. The program may be used for a Tetrahedron, 
Octahedron, or Icosahedron, depending upon the coordinates 
chosen as input data. The output is given in units based 
upon a radius of 1 for the spherical form and therefore, 
may be used as a basis for determining large structures. 

The example of input data is given in Table 1.6. The 
example of output data is given for a six frequency Icosa- 
hedral sphere, and may be read as in Table 1.7. The output 
takes advantage of symmetries within the spherical Icosa- 
hedron as discussed in the text material. Figure 1.56 
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THE COMPUTER PROGRAM DESCRIBED ON 
PAGES 1-97 to 1-144 
IS AVAILABLE FROM COSMIC 




1 .8 METHODS 4 & 5 


This mathematical and computer model was written for 
subdivision of a tetrahedron, octahedron, or icosahedron 
circumscribed by a unit sphere. The Icosahedron was chosen 
as an example to illustrate the geometry of the model. The 
polyhedron is oriented in a three dimensional rectangular 
coordinate system so that the vertices of one PPT are: 



(X2 .y 


(0, .850651, .525731) 



^ ^ 3 3 * ^ 3 ^ 


( .525731 , 0 , .850651 ) 
, 1 0 
"yryT 


- ( .850651 , .525731 , 0 ) 

where : t = 1 + /T 

2 
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with the intersection of the axis x, y, z located at the 
origin (0,0,0) of the polyhedron. Figure 1.58. 






Figure 1.58 

Due to the methods in which the three-way grids are generated 
only even frequency subdivisions may be used. 

Subdivision and Gridding for Method 4 
The PPT is subdivided into equal parts so that the 


vertices take the form. 


(xi + I 


- Xn , 


T yi + I 


Y2 - Vl > ^2 ■ \ 

^ ^ 1 ^ N / 


\ 

where N is the frequency of the structure and I is an interger 
such that 0 < I .< N. The values of I are unique for each vertex 
,nd are used to identify each vertex for each side as shown in 


F i gure 1.59. 


I-lh-6 







$2 



X,,Y, , Z, 

Figure 1.59 


Having found the unit division along the principal side of 
PPT, the points of subdivision are connected so that a three- 
way grid is generated with the lines of the grid perpendicular 
to their respective sides. Figure 1.60 



Fi gure 1 . 60 


I-lh-T 







Solve the equations for the x & y coordinates of the inter- 
section of PiP2 with PTPj^ 



D 


find the z coordinates 

For PjP^ let; (z^-z^ ) = 

(Yi- Y2) = 

^1 (Yi-Yz) + Yi (z 2 - Z3 ) = c^ 

For P3P4 let: - Z3) = a^ 

(Y 3 - Y4) = t>2 

- y^) + y 3 (z^ - Zj) = 

The intersection of PTF^with F7F'4& FTP^ are coincident and 
need not be found. 

All other points of intersection of the three-way grid are 
found in like manner and arestoped for final translation to 
the surface of the circumscribed sphere. 


1 .28 


1 .29 
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Subdivision and Gridding for Method 5 
With the following equations the planes consisting of 
the edge of the PPT and the origin (0,0.0) are rotated from 

3-space , Figure 1.61. 


Z 



X" = AjX + yiY + ViZ 

y = X 2 X + 

z = X3X + UgY + V3Z 


Where x. 


Z " - a X i 3 


X 


2 ’ 


y, V are direction cosines of the X"-axis 

respectively wit h respect t o the old axis 
/ ~ 2 2 

Xj ' /Xj/ Xi + yi + 7-1 

f ^ 2 2 

Ml =/yi/ ^ 

/ ^ 2 ~”^ z ^ 

V'l = yzi/ Xj + Yi + Zj 

va- V 3 similarly. 


, Y^-axis, and 
and are found by 
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Due to the method in which the three-way grid is generated 
only even frequency subdivisions may be used. The principal 
side of the PPT is subdivided into equal arc units by the 
following method: Figure 1.62 and 1.63. 

FIND: the angle contained within the rotated triangle 

consisting of PiP2» and the origin with the vertex located 
at the origin. 

cp = Arctan ( ^y ^ ) • r 1 .31 

— 

X 

Where r = 1 and is considered constant 


Y 



Figure 1.62 

THEN: subdivide the angle 4, into N angles 

® ^ • T 1.32 

N 

where T = Increment 1 to N 
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Figure 1.63 


The points of intersection of OP3 and found; 

PTP2 is y-Yi = ^2"^! 

X - X j X2 - Xi 

0^3 isy-0 =y3"0 

X - 0 X3 0 

The equation takes the following form: 

^2 is x(yj, - Yi) + y(Xi - Xj) = yi(xi - X2) + X (Y2 

let (Y 2 " y 1 ) = 

(X2 “ Xj) - 

y,(xi - X 2 ) + Xi(y 2 - Yi) = Cl 



Solve the equations for the point of intersection: 



1 .34 


Rotate the points of intersection along the PPT edge from 
2-spaces back to 3-spaces. 

X = AjX " + y + v^z 1.35 

y = A^x" + y^y" + V2Z 

Z = A^x " + ygY " + V^Z 

where A, y, v, are direction cosines of the X"-axis, Y"-axis, 
and Z"-axis with respect to the old axis and are found: 



A , A ; yp, y ; and V2, v„ are found similarly. 

Retain the co-ordinates along the edges and as shown 

in Figure 1.64. 
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Xi,Y,,Z, 

Figure 1.64 


After finding the unit divisions along the principal sides of 
the PPT, the points of subdivison are connected thereby creating 
a grid network. Since the units along the principal sides are 
not of equal length, the gridding will create "windows". The 
centers of these "windows" must be found to establish the 
three-way grid on the PPT. Figure 1.65. 




Figure 1.65 

The gridding and windows are found by the following 
method: From the coordinates along the ed'ges of S^, S2 and S3, 

calculate coordinates of the window by finding the intersection 
of F7P2 with P7F4 and P7P2 with and with by using 

the two point form of the equation of a line in three-space for 
the three lines and solve simultaneously for the points of 
i ntersecti on . 
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Find the z coordinate: 


For P ;^P 2 let: 


For P P let 

3 4 


(Za-Zj) = 

(yj-y^) “ 

z(y-y)+y(z-z)=c 
(z -z ) = a 

^43 2 

(y -y . ) = 


S(y3-y4> ^ y3("4-"3> = s 

The other two vertices of the window are found in a similar 
manner. Once the coordinates for the vertices of the window 
are determined, its center is found by the following method: 
cw=(w+w+w)/3 

1 2 o 

where: cw = center of the windows 

w = the X, y, or z coordinate of the 
vertices of the window. 



Translation of the Grid for Method 4 & 5 

The vertices of the three-way grid are then translated 

to the surface of the sphere along a line passing through the 

respective vertex and the origin (0,0,0) of the polyhedron by: 

X = rx 
1 1 

d 

y/ = ry^ 

~ 1 .37 

z = rz 
1 1 

“d“ 
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Where : 
and 
and 
where : 


/2 2 2 

'xi + yi + 

distance from origin to P 


= the radius of the unit sphere 
= 1 


Using the translated coordinates, this program finds the 
lengths of the elements of the structure {a), the angle between 
pairs of elements (face angle a), the angle between the elements 
and a radius from the origin to an endpoint of the element 
(axial angle fi), and the angle between adjacent faces of the 
structure (dihedral angle g). Figure 1.66 
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Figure 1.66 

To find the angle between elements the face ^ a, we use the 

coordinates of their endpoints, Thfe vertex of the angle is a 

common endpoint to each element and is translated to the origin. 

The other two endpoints and are translated in the same 

manner. Letting (x,, , zj and (x,, , z^) be the points 

resulting from the translations of the endpoints P and P 

1 2 ’ 
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cos a 


S ><2 ^ ^ ^"2 


d d 
1 2 


where d^-yxi + Yi 

and = \/^2 * ^ 2 * ^2 


1 .38 


a is the desired angle. 

To find axial angles the above method is used except that 
the vertex is established at one end of an element and the 
origin is used with the other endpoint to define the angle. 

The desired angle is Q.. 

The angle between two adjacent faces, the dihedral 3, 


is found using 


cos 3 = 

^1^2 ^ ^1^2 ^1^2! 


V 

/aj^ + + Cj^/a 

2 2 ^2 
+ B + c 
222 


1 .39 


where 


3 is the desired angle. 

A X + B Y + CjZ + Di = 0 defines the plane containing one 
face aid B,Y . C,Z + D, = 0 defines the plane containing 

the other face. The negative sign is used because the obtuse 
angle is desired. 

The A, B, and C for each plane are computed as 


A = 



B 


X 

1 

X 

2 

X 

3 


z 

3 


z 

2 

z 

3 


1 

1 

1 

1 

1 

1 
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1 


c = 


^3 Y3 1 


Where (X^, , Z^), (X2, Y2 . 1 ^) , and (X^, , Z^) lie in 
the plane. In particular the three vertices of each face 
are used. 

The length of the elements are found by using the 
general equation: 


£ = 


2 2 

(P - P ) + (P - p ) + (p - p ) 


1 


V 

" 2 


1 .41 


£ is the desired length. 

To reduce total output, this program takes into account 
certain symmetries and outputs only a part of the total angles 
and lengths. The rest of the values are the same as at 
least one outputed value and can easily be found using the 
following symmetries. Figure 1.67, 
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THE COMPUTER PROGRAM DESCRIBED ON 
PAGES 1-163 to 1-184 
IS AVAILABLE FROM COSMIC 




1 . 9 Methods 6 & 7 


This mathematical and computer model was written for 
subdivision of a tetrahedron, octahedron, or icosahedron 
circumscribed by a unit sphere. The icosahedron was chosen 
as an example to illustrate the geometry of the model . 

The polyhedron is oriented in a three dimensional rec- 
tangular coordinate system so that the vertices of one 
PPT are: 




(X2 , y2 , Zz) 


(^3* ^3* Z 3 ) 


0 , j/r 

yr 


1 


/T/7. 

( 0 , .850651 , .525731 ) 

1 > 0 

yr / 5/7 

« ( .850651 , .525731 , 0 ) 

J » 0 , 

v/T/T /T 

- (.525731 , 0 ,. 850651 ) 


where: t = 1 + 

2 
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The intersections of the axis X, Y, T is located at 
the origin (0,0,0) of the polyhedron. Figure 1.68. 



Figure 1.68 
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The PPT is first translated by the following equations 
such that the origin lies at the first point 



Arctanl y 



Cos(T )Cos(T )+y Sin{T )Cos(T )+z Sin(T 



vision, only triangle ABC is subdivided. The remaining 
sections of the PPT are found through rotations and reflec- 
tions of the basic unit. 



Note: ABC is a 



Figure 1.69 

Subdivision for Method 6 

The Line AB is subd ivied into parts chosen as equal 
are divisions of the central angle of the polyhedron by the 
following equation. Figure 1.70. 



where: 9 = ^ / 2 6 

6 = the central angle of the polyhedron 
0 < I<N 

(xi, yi) and (x 2 > y 2 ^ represent any two 
points on the PPT edge 


0 


Figure 1.70 


Subdivision for Method 7 


The Line AC is subdivided into parts chose as equal 
arc divisions of an angle made up of the triangle AC and 
the origin of the polyhedron with the origin (0, 0, 0) being 
the center of the triangle of subdivision. The following 
equations are used for this subdivision. Figure 1.71 

Arcs in / ^ J (x^-X 2 )^ "^ ( y i ~ y 2 ^ ^ ^ 

2 (x^-x^) 

X = — — Sin / I A \ 

3y(X2-Xi)2 +(y2-Yi)2 N ) 

2 (y^-Yi) 

V = , Sin/ ^ ^ 

3y(x,-xJ2 +(y,-y,)2 V N } 


where: A- the angle between 

and ^ 

0 < I<N 

(xi,yi) and (X 2 ,y 2 ) repre- 
sent any two points on the 
PPT face bisector. 


(X2.Y2) 



Figure 1.71 

Gridding & Projection for Methods 6 & 7 

The points of subdivision for methods 6 and 7 are 
stored in a matrix and are used in the gridding process for 
each respective method. Figure 1.72 

Matrix PT(A,B,C) 

spheres A = the type of point 

1-on the edge of triangle 
2~on bisector of angle 

3- on internal points of triangle 

4- external points 

B = line value and x, y, z values 

1-3 represents x, y, z values of line 1 

4-6 represents x, y, z values of line 2 

7-9 represents x, y, z values of line 3 

C = the number of points on the line. 


I-I9I 



\c\ 

2C1 



Figure 1.72 


Once the subdivisions are found, they are used to find 
points of divisions on the other two sides of the right 
triangle by the following equations. Figures 1.73 & 1.74. 

X = (Yj +Xj [MjiD - ( V3 -X3 [MjJ) 


y " x(M 2 i) + yi“Xi(M 2 ].) 

whereiM^^ =-X 2 -Xi (neg. reciprocal of slope) 

y2-yi 

^34 = ^ 3"^4 (slope) 



Figure 1.73 






Reflection of points in twor-space are found by: 

Figure 1.75. 

\ = [(yj - (yj 

(Mj -Mj) 

y^ " ^ 4 '^! '^3 

where 1 • 50 

M2 = (y^ -y,) 

(X 2 -Xj) 

and the reflected points are: 

X = Xj + 2 (x^ - Xj ) 
y = yi + 2 (y^ - yj 


i (Xj.Yj) 




(X.Y) 


(X 4 ,Y 4 ) 


4(X3.Y3) 


Fi gure 1.75 

Rotations of points in two-space are found by: 
Figure 1.76. 

X = (x2"Xj) Cos $ - (y 2 -yi) Sin $ +Xj 
y = (X 2 -Xj) Sin $ + (y^-yj Cos $ +y^ 
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1 .51 



in matrix form 


(x,y) 

Cos 

^ - Sin $ 


Si n 

$ Cos $ 



Fi gure 1.76 


The points for the internal gridding are found by 
the equations: Figure 1.77. 


X = (BjCi - BjC^) 
(AjBj - AjBj) 


1 .52 
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where : 


y = (A^Cj - AjCj) 

"(AiBj - BjAj) 

*1= - >1 
Bi= Xj - 

Ci= Xj (yi - Ya) + Yj (x^ - x^ ) 
*2= ^4 - >3 

><3 - S 

^2"" ’'3 ■ >4) * >'3 (^<4 - X3) 




Fi gure 1.77 


Through rotations and reflections (equations 1,50, 1.51) 
of the basic unit, the entire PPT three-way grid is found. 
Figure 1.78 
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A 


B 


Figure 1.78 

The external points of the PPT are found by equation 
1.52 and are rotated into their respective plane by equation 
1.44 where: the angles of rotation are found by equation 

1.43. The external points are rotated by equation 1.51 

to their respective positions. 

All points of the three-way grid are then rotated into 

three-space using equation 1.45 
where: = 2n - 

T = 2n - T 

y ^ 

T = 2n - T 

X X 

are the angles of rotation. 

The origin is then retranslated to its original 
position by: 


.53 
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X 


= X + T 


= y + T, 
= z + T. 


1 .54 


where: (T^, , T^) 

the coordinates to which the origin was 
originally translated to. 

All points of the three-rway grid are then projected 
to the surface of the sphere by: Figure 1.79 


/ 

X 

/ 

y 

z' 


X 

Di s 

y 

Di s 
z 

Di s 


1 .54 


where: Di s = v/(x) ^ + (y) ^ + (z) ^ 



Fi gure 1.79 
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For the PPT, the number of: 
edges = 3v(3v - 2) 

8 

half edges = 3v 
2 

faces = 3v(v - 2) 

half faces = 3v 
vertices = 3v(v + 2) 

+ 1 

8 

where:v = frequency and must be even 

For the total spherical form, the number of: 

r ^ 

Edges = 3 ^ 

2 

Faces = 

~ 2 ~ 

2 

Vertices = + 2 

4 

where: E = no. of edges in polyhedral unit. 

Using the coordinates, the lengths of the elements of 
the structure (z), the angle between pairs of elements 
(face angle a), the angle between the elements and a radius 
from the origin to an endpoint of the element (axial angle q) , 
and the angle between adjacent faces of the structure 
(dihedral angle 3 ), are calculated. Figure 1.80 
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Figure 1 . 80 


To find the angle between elements^ the face aj. a, we use 
the coordinates of their endpoints. The vertex of the angle is 
a common endpoint to each element and is translated to the 
origin. The other two endpoints, and are translated 
in the same manner. Letting (x^, y^, Zj) and y z^) 

be the points resulting from the translations of the end- 
points P^ and P^ , 


cos a 


where d^ 


XjX^ + + 21^2 

d id 2 
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and 


= v/x2^ + 


a is the desired angle. 

To find axial angles, the above method is used except 
that the vertex is established at one end of an element and 
the origin is used with the other endpoint to define the 

angle. The desired angle is 

The angle between two adjacent faces, the di hedral 3 » 

is found using 

® , , 1.58 

7Aj2+ Bj2+ y Bj2+ 0^2 

where 

3 is the desired angle. 

A^X + BjY + CiZ + Di = 0 defines the plane containing 
one face and 

X + BoY + CoZ + D. = 0 defines the plane containing 
the other face. 

The negative sign is used because the obtuse angle is desired. 
The A, B, and C for each plane are computed as 


A 


Y 


1 



1 

1 

1 


B 


X3 


1 .59 
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c = 


where (X Y Z ). (x Y 2 ). and (X . Y . Z ) lie in 

222 333 

the plane. In particular, the three vertices of each face 
are used. 

The length of the elements £ are found by using the 
general equation: 

£ = 


^1 ^2 
^ is the desired length 


- P )2 + (P _ p )2 + (p - p ) 


yi 


y2 


At the time of publication, the computer programs for 
these methods had not been completed and therefore, have not 
been included in this report. 
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